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Abstract 

In this Reply, we respond to the above Comment. Our computation 
[Phys. Rev. D 91 (2015) 074512] only took into account pure QCD effects, 
arising from quark mass differences, so it is not surprising that there are 
discrepancies in isospin splittings and in the Sigma - Lambda mixing angle. 
We expect that these discrepancies will be smaller in a full calculation 
incorporating QED effects. 
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We agree with Gal [T] that lattice numbers for the mixing angle are not 
complete until numbers for the isospin splittings are correctly accounted for. In 
our paper, [2] the calculation only includes QCD effects, arising from quark mass 
differences. For a complete result, a full calculation with QED effects added is 
also needed. 

Before addressing the main issue, we would hrst like to mention that it is 
not surprising that our results satisfy the Coleman-Glashow [3] and Dalitz-von- 
Hippel (eq. (3) of [1]) relations. Our £t function has all the SU{3) symmetry 
constraints built in, so it automatically obeys every symmetry relation to the 
appropriate order. In this case, violations of the Goleman-Glashow relation are 
at 0{6m^), while violations of the Dalitz-von-Hippel relation are at 0{6m‘^). 

We believe that lattice QGD gives reliable numbers for the part of the isospin 
splitting due to quark mass differences, and that the differences between the 
values in [2] and experiment are mainly due to electromagnetic effects. Our 
reasons to believe that lattice numbers for the purely QGD contribution are 
accurate are: firstly the lattice gives good values for the splittings between the 
multiplets (W, A, S,S), and it is hard to see how there could be a systematic 
error that would spoil the isospin splittings without also showing up in splittings 
between multiplets. Secondly at leading order, pure QGD relates the isospin 
splittings to the splittings between the multiplets through the relations 
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(see Fig. 2 of |5]). Our simulations show that this leading order formula only 
has minor corrections from higher order effects, and the above relations hold 
reasonably well, with a value ~ 0.022 for the quark mass ratio (m^ — rriy)/{mg — 
rriy)- The main systematic uncertainty in this mass ratio are due to the difficulty 
of correcting for electromagnetic effects in the pseudoscalar meson sector. In 
Fig. [1] we show the splitting values. The squares and circles are consistent with 
each other, but do not reproduce the full result without adding a non-QGD force 
(namely QED). 

The relations ([T]) do not hold in the real world, we take this as evidence that 
the QED corrections to the isospin splittings are substantial. Simulations with 
QED included |6] show that QED effects on the isospin splittings are comparable 
with the effects of the (m^ — m„) difference, and that the combined simulation 
reproduces the physical numbers very closely. Unfortunately we do not have 
lattice results for the QED shift in the mass, so we can not estimate the QED 
effect in the mixing matrix element. 
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Figure 1: Baryon octet (mass splittings)^. Black squares are from eq. ([T]), while blue 
circles the lattice QCD results. Both use the pure QCD numbers from table 5, [2]. Red 
stars use the physical values (i.e. for QCD and QED). 


In the paper by Isgur [7] there are electrostatic (Coulomb) contributions to 
the isospin splittings, but they cancel completely for the combination [(M^o — 
^E+) ~ {Mn — Mp)], which, by the Dalitz-von-Hippel relation, is proportional to 
the S-A mixing angle. If this holds, we would expect to see some shifts in the 
isospin splittings, but no Coulomb contribution to the mixing angle. However, 
this exact cancellation is model-dependent. The Coulomb contributions in [7] 
are calculated for a particular Gaussian wavefunction, which has more symmetry 
than required by QCD. For equal-mass quarks the Isgur spatial wavefunction is 
completely symmetric under exchange of any quark pair. Since the octet is a 
mixed symmetry multiplet, we do not expect complete symmetry. In the proton 
{uud) there is no theorem that says (l/r„„) = {l/r„rf). 

We have applied our flavour analysis to QED effects too. We hud that there 
are five allowed coefficients for the electromagnetic effects, one of which just shifts 
all masses by the same amount, so making no contribution to splitting or mixing. 
In terms of these coefficients, we hud (assuming small mixing angle) 


= 
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Mlo = ^[8B^+ 8Bf^+ 2B™+ 9Bf^+ 3B™]e^ 

M|+ = -[2B™+ 3B^ - B™+ 3Bf^]e^ 
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M|_ = -[2B^ + Bf^]e^ 
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M|o = ^[2Bf^+ 2B™ + Bf^+ 3B™]e^ 

as the general expression for the electromagnetic contribution to the masses, and 

(S0|M^|A) = ^ [-2i?r - 351“ + 35r] (3) 

for the mixing matrix element. (We £t to the squares of the masses - the form 
of the expansion is of course the same if a £t to the masses themselves is made.) 
It may readily be checked that these electromagnetic contributions satisfy the 
Coleman-Glashow and Dalitz-von-Hippel relations, for any value of the B™ co¬ 
efficients. 

In the Isgur model, these coefficients are not all independent. The Coulomb 
terms of the Isgur model follow the pattern 

= Bf^ = -Bf^ , B™ = 0 (4) 

which in turn ensures that there is no Coulomb contribution to the mixing, see 
eq. ([3])- However, these inter-relations are model dependent, and we need lattice 
calculations to see how well they hold. 

We are currently working on a combined simulation with QED effects in¬ 
cluded. Preliminary results suggest that QED effects do account for much of the 
difference between the current QCD-only results and the experimental values. 
In the Isgur model it is expected that although the individual splittings will be 
changed by QED effects, the QED contribution to the E-A mixing angle will 
cancel. This is however a model dependent statement, with other choices for the 
spatial wavefunction the cancellation would not be complete. Our joint QED and 
QCD results are not yet at the point where we can comment on how much the 
mixing angle is changed by QED, but we are grateful to Gal for bringing the issue 
to our attention, and will certainly include a discussion of the question when we 
have our hnal results. 
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